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ABSTRACT 

Micro particles are today a major technological and economic issue. They enable very promising innovations in 
many fundamental areas such as: health, energy or industry. Despite the numerous studies devoted to hydrodynamics in 
agitated tanks, there is still a glaring lack of information on the temporal evolution of fluid flow properties. The 
knowledge of this dependence has been often crucial , for example, in the food , chemical , biochemical or petrochemical 
industry is carried out in tanks or in mechanically agitated reactors. The micro particles are colloidal solutions composed 
of particles of micrometric size suspended in a liquid. Their astonishing thermal properties have been the subject of 
intense investigation during the last decade. It is proposed in this work to study the influence of the concentration of 
micro particles (al2o3, cuo, etc.) Dispersed in a base fluid on the hydrodynamic fields of fluid in flat-bottomed cylindrical 
and non-quenched fluid provided with a stirrer. Anchor type was undertaken by numerical simulation using the fluent 
cfd 6.2.13 based on the finite volume discretization method of the navier stokes equations (uvp) equations and the energy 
equation. We aim at a synthesis of the simulated system by comparing the power consumed by the micro particle 
technique. 
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1. INTRODUCTION 

In recent years a new type of fluid is emerging: the nanofluid. The nanofluid is a coolant such as water, 
the oil to which is added nanoparticles of metal (Al, Cu, Ag, Au, etc.) or metal or non-metallic oxide (Si0 2 , A1 2 0 3 , 
Ti0 2 ), and other (allotropic forms of carbon) in relatively small amounts (0.1 to 6% of the total volume) [1]. 

Nanofluids are a new fluid category that was introduced for the first time by Choi [2, 3]. Indeed, it 
involves dispersing nanoscale particles in conventional fluids commonly used, in order to improve the properties of 
the latter. Nanofluids are generally classified according to the nature of the nanoparticles that compose them in 
three categories, namely, metal, metallic and non-metallic oxide [4, 5]. Nanofluids find their applications in several 
sectors, including domestic, engineering and biomedicine. 

The development of nanotechnologies is a challenge for the years to come [6]. Nanoscience studies the 
structure of matter at the molecular level, therefore at the nanoscale, with the aim of innovating in many fields such 
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as: the environment, the production and storage of energy, the cooling of electronic components, biology, medical 
diagnosis, water treatment, etc [7, 8, and 9]. 

Microparticles will be today a major technological and economic issue. They enable very promising innovations 
in many fundamental areas such as: health, energy or industry [10, 11, and 12]. Despite the many studies devoted to 
hydrodynamics in agitated reservoirs, there is still a glaring lack of information on the temporal evolution of the flow 
properties of fluids [13, 14, and 15]. The knowledge of this dependence is often crucial for example in the food, chemical, 
biochemical or petrochemical industry is carried out in tanks or in mechanically agitated reactors [16, 17, 18, and 19]. The 
micro particles are colloidal solutions composed of particles of micrometric size suspended in a liquid. Their amazing 
thermal properties have been the subject of intense research over the last decade [20, 21, 22, and 23]. 

The main objective of the present work is to study the influence of the concentration of micro particles (A1 2 0 3 , 
CuO, etc.) dispersed in a base fluid on the hydrodynamic fields of the fluid in a cylindrical and non-hardened flat-bottom 
fluid provided by an agitator. The type of anchoring was done by numerical simulation using the CFD 6.2.13 method based 
on the finite volume discretization method of the Navier Stokes equations (UVP) and the energy equation. We aim at a 
synthesis of the simulated system by comparing the power consumed by the microparticle technique. 

Many industries (chemical, pharmaceutical and food) use, produce or transform complex fluids. These are high 
viscosity fluids, often with non-Newtonian behaviour. Their preparation generally requires a heating step followed by a 
cooling step. Heating is essential to facilitate manufacture (chemical reaction and / or physical-chemical transformation) 
[24, 25, and 26]. 

The cooling or the finishing step consists in bringing the fluids to a temperature close to ambient temperature in 
order to approach the conditioning step. The cooling is most often carried out in the same stirred tank [27]. 

There are two strong demands in the industry: firstly, ensure the quality of finished products. Secondly, work in a 
more and more concentrated environment for pollution and cost problems, in order to know and check all the involved 
aspects in the process, in particular hydrodynamics and heat transfer [28, 29]. 

2. NUMERICAL MODELS 




Figure 1: The Geometry of Mixing System 

The mixing systems are presented in figure 1. The tank is a flat-bottomed vessel (inside diameter: D) equipped 
with an anchor agitator (diameter d). H is the fluid height. Values for the agitator width ‘L’, agitator-to-wall clearance ‘w\ 
impairs thickness ‘e’ and geometrical characteristics are reported in Table 1. 
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Tablel: Geometrical Characteristics of the Agitator 


D 

d 

da 

e 

L 

w 

1 

0.96 

0.023 

0.027 

0.067 

0.04 


3. MODEL EQUATIONS AND COMPUTATIONAL PROCEDURE 

The equations describing fluid flow are derived from the conservation of mass and momentum which are also 
known as Navier-Stokes equations for fluids of constant physical properties [1, 30]: 

Continuity: 


divV = 0 


( 1 ) 
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Density: 

Pnf = <PP S + (l - P)p f 

And 

The viscosity of the nanofluid is given by Brinkman [1]: 


Pnf = 


Pf 


The variables of the previous equations are dimensional as follows: 
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In dimensionless form: 
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The Reynolds number is defined by 
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Finally the equation of the amount of movement along the X and Y axis is: 
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U\ 
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3.1. Boundary Conditions 

In the CFD code used two reference frames can be treated. One is stationary, the impeller is in rotation and the 
tank is fixed. The other frame, which is rotating, the impeller is kept stationary while the outer wall of the vessel is given 
an angular velocity equal and opposite to the velocity of the rotating frame. Constant boundary conditions have been set 
respecting a fixed reference frame (RRF) approach. The boundary conditions for velocity are fixed on the gate impeller and 
the vessel. On the gate impeller: Moving reference frame implemented in the CFD code. 

3.2. Validation 


To verify the model of the fluid and to validate our results we compared our power number with previous works 
on an anchor agitator. Figure 2 gives an example of comparison between our results and a numerical work of Marouche 
[15]. The results show a satisfactory agreement. 



Figure 2: Variation of Power Number as Function of Reynolds Number 


4. RESULTS AND INTERPRETATIONS 




Figure 3: Tangential Velocity on the Blade and Figure 4: Tangential Velocity on the Mediator 
its Extension according to Re for C= 20% according to Re for C= 20%. 
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Figure 5: Radial Velocity on the Mediator 
According to Re for c= 20%. 


Figure 6: Radial Velocity on the Blade and its 
Extension According to Re for C= 20%. 




Figure 7: Tangential Velocity on the Blade and 
its Extension according to Concentration C 


Figure 8: Tangential Velocity on the Mediator 
according to Concentration C 




Figure 9: Radial Velocity on the Blade and its Figure 10: Radial Velocity on the Mediator 
Extension as a Function of Concentration C as a Function of Concentration C (CUO) 
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Figure 11: Tangential Velocity on the Blade and its Figure 12: Tangential Velocity on the Pale Mediator 
Extension as a Function of Concentration for and its Extension as a Function of Concentration for 
Re = 415 (A1 2 0 3 ) Re = 415 
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Figure 13: Variation of Np as a Function of Re 
with a Concentration of 20% for (A1 2 0 3 ) 


Figure 14: Variation of Np*Re as a Function of Re 
with a Concentration of 20% for (A1 2 0 3 ) 




Figure 15: Variation of Np according to Re Figure 16: Variation of Np * Re according to 
for different Concentration (A1 2 0 3 ) Re for different Concentration (A1 2 0 3 ) 
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Figure 17: Variation of the Tangential Velocity Figure 18: Radial Velocity Variation on the Blade 
on the blade and its extension for (A1 2 0 3 ) and its extension for (A1 2 0 3 ) 

and (CuO) for Re=415 and (CuO)for Re=415 



Figure 19: Variation of Np as a Function of Figure 20: Variation of Np * Re as a Function 
Re for (A1 2 0 3 ) and (CuO) of Re for (A1 2 0 3 ) and (CuO) 


Les Figures 3, 4, 5, 6 and 7 show the variation of the tangential and radial velocity of the agitator and its extension 
and on the mediator according to the radial position. We studied the influence of Reynolds number on the nature of the 
flow by varying this parameter from 0.7 to 415 for nanoparticle concentration C=20%. The tangential component of the 
velocity varies in the following manner (Figure 3). 


The tangential component of the velocity varies as follows: the growth of the tangential velocity of the shaft of the 
agitator at the end of the blade a linear increase between its two ends, once completed, a sudden decrease in the velocity of 
the wall of the tank which is null. It is noted again that the tangential velocity increases with the increase of the Reynolds 
number. 


It is noted that upstream of the blade the tangential velocity Vt* is negative for the low Reynolds numbers, when 
Re is large enough, the tangential velocity is always positive, against the radial velocity Vr*: this means that the end of the 
blade sucks the fluid from the center of the wall of the tank. On the mediator Figure 4, the velocity is close to zero for low 
Reynolds numbers, but reaches a maximum with increasing relative inertia. Moreover, this maximum is lower than the end 
of the blade, which means that the flow rate in this zone is delayed beyond this maximum speed decreases as one 
approaches the tank, although more inclined than that from the ascending part of the curve. 
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The influence of the Reynolds number on the profile of the radial velocity Vr* on the blade Figure 5, we see 
thatthe variation of the radial velocity is perceptible Clearly between the axis of the agitator and the blade, this variation is 
rapidly decreasing at a low Reynolds number, we noted the existence of negative velocities which lead to the formation of 
the zone Recirculation at the vicinity of the blade and near the tank where the flow is delayed in this area. In addition, on 
the mediator Figure 6, the radial velocity is zero. 

Figures 7 8 9 10 show the variation of the tangential velocity and its extension and on the mediator for different 
concentrations for a flow regime Re = 415 and a nanofluid (CuO), we note that there are no important variations because 
these quantities are local, similarly for the fluid nano (A1 2 0 3 ) figure. 11, 12. 

We note that the power number varies linearly with the inertia (Re) on a logarithmic scale, produced by the 
against Np*Re, the laminar flow constant is Re <276 equal to 800 for our mobile, all above this value, there is Np * Re, 
which means that energy consumption Reynold's number (figure 13, 14). 

Figures 15 and 16 illustrate the variation of the power number Np as a function of the inertia and the product 
Np*Re as a function of the concentration for a nano fluid (A1 2 0 3 ), there is a slight variation for low Reynolds numbers. As 
soon as the inertia is increased this variation becomes important, this can be explained by the large viscous forces in the 
case of nanofluids and is due to the improvement of the thermal conductivity of the nanofluid, when the volume fraction of 
the nanoparticles increases. 

There is no variation of the tangential velocity on the blade and its extension and on the mediator between the two 
nanoparticles (A1 2 0 3 and Cuo) figure 17, 18. But the variation is remarkable in terms of number of power. 

Figures 19, 20 represents the variation of the Reynolds number as a function of the volume fraction of the 
nanoparticles for different nanofluids and the product Np* Re. Note that the number of power Np for CuO is greater than 
A1 2 0 3 ; this variation is remarkable for large Reynolds numbers. 

5. CONCLUSIONS 

The hydrodynamics of nanofluids in complex geometries such as stirring systems are modestly discussed. 
However, it plays an important role in the control and the effectiveness of the treatment because it conditions the transfer 
of heat and mass. 

The problem is analyzed numerically. The latter is not treated by a code, but by using an existing industrial code 
which is the current code. The numerical simulations carried out in 2D geometries allowed to show significance by the 
rheological behavior of the fluid. 

We used two types of nanofluids (A1 2 0 3 + water), (CuO + water), each with different thermo physical properties. 
We have assumed that nanofluids can be assimilated to homogeneous fluids. 

From our study, we have led to results that differ significantly from the results found by the many authors 
mentioned above. 

It is true that the addition of (A1 2 0 3 , CuO) nanoparticles, in pure water increases its thermal conductivity and thus 
improves the conduction heat transfer in the nanofluids and also increases the energy consumption which is dissipated from 
the agitator to the fluid. 
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The number of power is higher in the nanofluids compared to the base fluid and it has also been found that this 
number in the nanofluids increases with the increase in the concentration of the nanoparticles. 
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NOMENCLATURE 


d 

: Impeller diameter, m 

D 

: Rate of strain tensor, dimensionless 

D 

: Vessel diameter, m 

H 

: Impeller height, m 

Hv 

: Vessel height, m 

da : 

: Shaft diameter, m 

e : 

Impeller thickness, m 

N : 

Rotational speed, rd/s 

Np : 

: Power number, dimensionless 

P : 

Power consumption, W 

Vr 

: Radial velocity, m/s 

Vt : 

Tangential velocity, m/s* 

W 

: Agitator to wall clearance, m 

d, reg 

: Regularization parameter of Bercovier andEngelman 

Bi 

: Bingham number, dimensionless 

He 

: Hedstrom number, dimensionless 

Re 

: Reynolds number, dimensionless 


Greek Symbols 


rjap 

: Apparent viscosity, Pa*s 

rjO 

: Cut off viscosity for low shear rate, Pa*s 

h 

: Limiting viscosity at infinite shear rate, Pa. s 

P 

: Density, kg/m3 

T 

: Stress tensor, Nm 2 

To 

: Yield stress, Nm 2 

7c 

: Shear rate, S-l 
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t : Cut off shear rate, Bingham model 
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